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Abstract. We study neutralino dark matter (DM) with large singlino fractions in the next- 
to-minimal supersymmetric Standard Model (NMSSM). We perform a detailed analysis of 
the parameter space regions of the model that give rise to such singlino-dominated neutrali- 
nos while satisfying the constraints from Higgs boson searches at the Large Electron Positron 
(LEP) collider and the Large Hadron Collider (LHC), as well as from 6-physics experiments. 
We find that this DM can yield a thermal relic density consistent with the Planck measure¬ 
ment in mass regions where the lightest neutralino of the minimal supersymmetric Standard 
Model (MSSM) generally cannot. This is particularly true for lighter DM masses, either less 
than 10 GeV or between 60 — 100 GeV, and for heavier DM masses, between 500 — 1000 GeV. 
We then analyse the prospects for indirect detection of such DM at the IceCube neutrino 
telescope, assuming the complete 86-string configuration including DeepCore. We also con¬ 
sider the added sensitivity to low-mass DM with the proposed PINGU extension. We find 
that IceCube is sensitive to some regions of the NMSSM parameter space containing singlino- 
dominated DM and that a subset of such model points are already ruled out by the IceCube 
one-year data. IceCube will also be sensitive to some parameter space regions that will not 
be probed by the upcoming ton-scale direct detection experiments. Moreover, we find that 
PINGU will be sensitive to DM in the 10 GeV range. 
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1 Introduction 

The nature of the dark matter of the Universe is one of the major unresolved problems in 
physics. Supersymmetric models with unbroken R-parity provide a possible candidate for a 
weakly interacting massive particle (WIMP), the lightest neutralino, which has been studied 
extensively (see e.g. [1-4]). There are good phenomenological and theoretical reasons apart 
from the DM problem to consider such supersymmetric (SUSY) extensions of the Standard 
Model (SM), including amelioration of the hierarchy problem or the fine-tuning of the Higgs 
boson mass and unification of the gauge couplings. 

The most commonly studied SUSY model is the MSSM [5, 6], which exhibits all of the 
above benefits. The price for these benefits is a common mass term (called the ;U-term) for the 
two Higgs doublets of the MSSM. This leads to another fine-tuning problem, namely why the 
^ parameter should lie at the electroweak (EW) scale, as is required by phenomenology, when 
its natural value would be the Planck scale. This is known as the //-problem. Furthermore, 
the mass of the lightest Higgs boson in the MSSM at the tree level has an upper limit equal to 
the Z boson mass. To fulfil the experimental constraints, large corrections to the Higgs mass 
from the stop (and top) loops are needed, which require a rather large stop mass, leading to 
an additional source of fine-tuning in the model. If we take fine-tuning arguments seriously, 
we may consider exploring alternatives to the MSSM. 

One such alternative is the NMSSM [7-10], which has been reviewed in Refs. [11, 12]. In 
the NMSSM, in addition to the two Higgs doublets of the MSSM there is a singlet Higgs field, 
which is the scalar component of a chiral singlet superfield added to the MSSM superpotential. 
The reason for introducing this additional scalar is that an effective //-term is now dynamically 
generated, so that the fine-tuned parameter // is no longer needed (and can be removed by 
imposing a discrete symmetry). This new scalar mixes with the other scalars from the two 
doublets, leading to a Higgs sector with seven Higgs bosons, compared to the five present in 
the MSSM. In addition, the fermion component of the singlet superfield, the singlino, mixes 


- 1 - 




with the higgsinos and gauginos, the fermion components of the Higgs doublet and gauge 
superfields, respectively, resulting in a total of five neutralinos. 

The interaction strengths of the neutralinos with other SM and SUSY particles are 
governed by their masses and compositions. In the MSSM, the neutralinos are broadly 
classified as gaugino-like, higgsino-like and mixed gaugino-higgsino states, depending on the 
relative sizes of the soft SUSY-breaking gaugino mass parameters, Mi^ 2 , and the fj. -parameter 
in the neutralino mass matrix. If Mi ^2 are much larger than /i, two of the neutralinos 
are predominantly higgsino-like and the other two are gaugino-dominated, while a smaller 
difference results in larger gaugino-higgsino mixing. Similarly, in the NMSSM, the fifth 
neutralino is either completely singlino-dominated or can be an admixture of the singlino 
and, chiefly, the higgsinos. This purely singlino or mixed higgsino-singlino state can be the 
lightest of the five neutralinos and hence the lightest supersymmetric particle (LSP) after 
diagonalisation of the mass matrix. If so, some interesting possibilities arise in the context 
of DM phenomenology, which are rather distinct from the MSSM. In particular, such a 
neutralino can give the correct DM relic abundance with a considerably smaller mass [13-17] 
than is generally possible for the lightest neutralino in the MSSM, given the LHC exclusion 
limits on SUSY. 

A light singlino-dominated neutralino LSP in the NMSSM has been the subject of 
several studies recently. It is an important candidate for the events near ~ 10 GeV [14, 15] 
reported by the CDMS H Collaboration [18]. With a mass near 35 GeV, such a neutralino 
can also explain [19-22] the ~ 1 — 4 GeV 7 -ray excess from the galactic centre reported by 
the Fermi Large Area Telescope (LAT) [23, 24] (for most recent results from Fermi-LAT, 
see [25]). Generally, the existence of a light DM candidate is also supported by some other 
direct detection experiments such as DAMA/LIBRA [26], CoGeNT [27, 28] and CRESST-H 
[29]. Experiments such as XENONIOO [30], LUX [31] or PICO [32] have released strong 
exclusion constraints on the interaction strength of a generic DM candidate with a nucleon 
for a considerably wide range of its mass. Signatures of a higgsino-singlino LSP that can be 
probed at the LHC have been studied in [33-37]. Indirect detection prospects for neutralino 
DM in the NMSSM have been studied in [38, 39]. 

In this article we revisit the neutralino sector of the NMSSM, focusing on the LSP 
with a significant singlino fraction, which can give sufficient relic abundance in the Universe. 
Through scans of the NMSSM parameter space, which were subject to the most important 
experimental constraints, we find its regions where such neutralino solutions are realisable. 
We then analyse the potential of the IceCube neutrino telescope [40] to discover or exclude 
such LSP solutions. The annihilation of neutralinos captured in the centre of the Sun results 
in a neutrino flux at Earth. These neutrinos interact with ice at the South Pole, producing 
leptons which emit Gherenkov radiation that can be detected by IceGube. We take into 
account the winter data from a complete 86 -string configuration of IceGube, i.e., including 
also the DeepGore strings, in order to perform a signal-to-background analysis of the good 
higgsino-singlino solutions obtained from our scans. The impact of the proposed extension 
PINGU [41] on the statistical significance for mainly low mass DM is also investigated in 
detail. 

The article is organised as follows. In section 2 we discuss the NMSSM and its neu¬ 
tralino sector. In section 3 we provide details of the model’s parameter space of interest, 
the procedure adopted for its scanning, the constraints imposed during these scans and the 
impact of these constraints. In section 4 we explain the neutrino flux from the Sun due to the 
annihilation of DM there. In section 5 we highlight some important features of the IceGube 
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neutrino observatory. In section 6 we present the details of the reconstruction of neutrinos 
from the Sun and the background for this signal at the IceCube. In section 7 we discuss 
the statistical treatment of the signal process and the sensitivity that can be reached at the 
IceCube after a given period of data accumulation. In section 8 we summarise our findings. 


2 The neutralino sector of the NMSSM 


The Za-symmetric CP-conserving NMSSM has no mass parameters in the superpotential, 
and is therefore also referred to as the scale-invariant NMSSM. It differs from the MSSM by 
the omission of the ^u-term from the superpotential and the introduction of a singlet chiral 
superfield S, while the Yukawa terms remain unchanged. The superpotential is written as 

IPnmssm = IPYukawa + XSHuHd + - S ^ , (2.1) 


where and are the two Higgs doublet superfields of the MSSM. The scalar potential of 
the NMSSM constitutes of the D-terms, which are the same as in the MSSM since S does not 
couple to the gauge bosons, the T-terms, plus the soft SUSY-breaking terms for the Higgs 
sector, which are given by 


Hsoft — 


m%_ 


I ff.. 


XAxSHuHd + -kA^S^ + h.c. 


( 2 . 2 ) 


The NMSSM scalar potential contains the new parameters X, k, A\, A,^, and ms compared 
to the MSSM. The parameters and m| are related to the Z boson mass, mz, 

tan/3 = Vujvdi where Vu and Vd are the vacuum expectation values (VEVs) of Hu and Hd, 
respectively, and s, the VEV of the singlet field, through the minimisation conditions of the 
scalar potential for EW symmetry breaking, and can therefore be eliminated. Here, instead 
of s we will use the effective /x-term, //eff = As, as a free parameter. The VEVs of the doublet 
Higgs fields also satisfy the relation v‘^ + v'^ = v'^ = 2rn^jg\ = (174 GeV)^, where mw is the 
mass of the W boson and 52 is the SU{2)l gauge coupling. 

The free parameters of the MSSM Higgs sector include m^, the mass of the CP-odd 
scalar, tan/3 and /r. In the NMSSM which is the 1x1 entry of the pseudo-scalar mass 
matrix before diagonalisation, gets traded for A\ and /x gets replaced by so that all in 
all we have six Higgs sector parameters (A, k, Ha, tan/3,/Xeff)- There are thus three more 
parameters in the NMSSM than in the MSSM, all of which originate in the Higgs sector. 

As 5 is a complex field, there are two additional physical Higgs bosons in the NMSSM 
compared to the MSSM. Eor a CP-conserving Higgs sector (as is assumed here), we have 
three CP-even neutral states and two CP-odd neutral states Hi and H 2 , where 

we take the states to be ordered in mass, so that Hi and Hi are the lightest scalar and 
pseudo-scalar states, respectively. 

The fermion component of S is called the singlino, S, which mixes with the gauginos, 
and IT 3 , and the higgsinos, H^ and H^. There are therefore five neutralinos in the 
NMSSM. At leading order the neutralino masses and mixings depend on the parameters of 
the neutralino mass matrix. If we introduce the vector = {—iB^,—iW^,H^,H^,S), the 
non-diagonal mass Lagrangian in the gauge eigenstate basis is given by 

^mass = M^O -h h.c., (2.3) 
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where Ai^o is the symmetric matrix 
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(2.4) 


with gi being the U{1)y gauge coupling. 

The neutralino masses and compositions at the tree level thus depend on the Higgs 
sector parameters X, k, fieg,Vu,Vd and the gaugino masses Mi and M 2 . The mass matrix in 
eq. (2.4) can be diagonalised by a unitary matrix N to give D = diag(m^o) = N*M.~oN\ 
for i = 1 — 5. If all the parameters in the mass matrix are real (which we assume to be the 
case here) then the matrix N is orthogonal, and we have D = NAt-oN'^ . The neutralino 
mass eigenstates are then given by Xi = 

The eigenvalues of Nij are all real, but can be positive or negative. (They can be 
made positive by a phase transformation.) They are not ordered in mass after performing 
the diagonalisation, but should then be reordered so that Xi is the lightest neutralino. The 
linear combination, 


^ = NuB^ + NuWl + + NuH^ + , (2.5) 


is thus our DM candidate. In order to describe the composition of Xi, we define the gaugino 
fraction in it as Zg = + |A^i 2 p) the higgsino fraction as = |iVi 3 p + |A^i 4 p and the 

singlino fraction as Zg = |iVi 5 p. As noted earlier, the focus of this study is a Xi with a 
non-negligible singlino fraction, which we define to be Zg > 0.05. Such a Xi could lead to 
some distinct phenomenological scenarios precluded in the MSSM. 

Let us note some properties of the tree-level neutralino mass matrix given in eq. (2.4), 
assuming the mass parameters Mi and M 2 of the gauginos to be very heavy so that they are 
decoupled from the 3x3 higgsino-singlino block. The 5x5 diagonal element, 2ks = 2Kfies/X, 
corresponds to the mass of the singlino. Thus, if this is small the Xi is more likely to be 
singlino-dominated after diagonalisation. Similarly, if is small, depending on the size of 
A, the lightest neutralino can instead have a larger higgsino component. To have a singlino- 
dominated WIMP, we therefore need a large /ieS and a small k. A large A also reduces the 
size of the 5x5 term, but at the same time it enhances the sizes of the off-diagonal terms, 
leading to a larger mixing. Note, however, that the presence of a certain amount of higgsino 
in the WIMP is necessary to obtain a realistic relic abundance. Thus should not be 
too large and A should not be too small. Evidently, the composition of such a WIMP is 
quite insensitive to the value of tan/3. In reality the situation is more complex, as the other 
elements of the matrix also need to be considered. Moreover, just as in the MSSM, this mass 
matrix is subject to radiative corrections [42-49] from various other sectors of the NMSSM. 
Thus, at higher orders the parameters of these sectors also become crucial. 


3 Model scans and constraints 

As noted above, beyond tree level the parameters of other SUSY sectors also need to be 
taken into account along with the Higgs sector ones, when drawing inferences about the 
neutralino sector. However, the most general EW-scale NMSSM contains more than 130 
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parameters in total. One way to reduce the number of free parameters is to consider a 
minimal supergravity-like scenario in which all the dimensional parameters are defined at 
the grand unification (GUT) scale and certain universality conditions are imposed on them. 
In this so-called constrained NMSSM (CNMSSM) (for some recent analyses, see, e.g., [50- 
52]) the scalar soft SUSY-breaking masses are unihed into a generic mass parameter tuq, the 
gaugino masses into and all the trilinear couplings into oq. Then, given m^, these three 

parameters along with the coupling A, taken as an input at the SUSY-breaking scale, MgusYj 
and the sign of /Teff constitute the only free parameters. All the parameters at MgusY are 
then obtained from these GUT-scale parameters using the renormalisation group equations. 

In order to allow more freedom in adjusting the Higgs and neutralino sector parameters, 
one can relax the above-mentioned universality conditions partially. In that case, similarly 
to the MSSM with non-universal Higgs masses (NUHM), mn^, and ms are disunihed 
from mo at the GUT scale and taken as free parameters, which can then be traded for k, 
/ieff and tan/3 at the EW scale. Also, a\ and a^, the GUT-scale input values, respectively, of 
the EW-scale parameters Ax and are disunihed from oq. This results in a total of nine 
input parameters and we refer to this model as the next-to-NUHM (NNUHM) here. 

Alternatively, without imposing the GUT-universality conditions, a convenient way to 
considerably reduce the number of parameters in the general EW-scale NMSSM is to assume 
that the matrices for the sfermion masses and for the trilinear scalar couplings are diagonal. 
Eurthermore, the soft SUSY-breaking parameters are taken to be real and those correspond¬ 
ing to the hrst two generations are unihed. Thus, with only about 25 or so parameters in 
total, one can study the most important low-energy characteristics of the model. In our 
analysis, while requiring the Xi to have Zg > 0.05, we also want to take into account all of its 
possible annihilation and co-annihilation channels resulting in the correct relic abundance. 
Some of these co-annihilation channels require the existence of a sfermion not much heavier 
than the y? itself. However, in order that the model candidate for the SM-like Higgs boson, 
Hobs) observed at the LHG [53-55] has a mass near 125 GeV, large corrections from the SUSY 
sector are necessary, implying sufficiently large sfermion mass parameters. Grucially though, 
these corrections are almost entirely dominated by the stops. In view of all these consid¬ 
erations, we impose the following universality conditions on the parameters of the general 
NMSSM: 


= = Mff = Mf, 

Q VI,2 Ol,2 2^1,2 

Mr = Mj = Mp , 

n ^1,2,3 ^1,2,3 


M2 = 

Aq = At = Ai, = Aj 


where Mq^ 2 , Mjj^ 2 soft masses of the squarks from the first two generations, 

and Me ^23 slepton masses, Mi^ 2,3 the soft gaugino masses and At^b^r the 

trilinear Higgs-sfermion couplings.^ This leaves us with a total of 14 free parameters to scan 
over. In the following we will refer to this model as the NMSSM-14. 

In order to obtain maximum possible configurations of the free parameters that yield the 
desired neutralino composition while satisfying important constraints coming from various 

^Note in particular that a unique soft mass has been used for all three generations of the sleptons. This 
is not the case for the squarks, for which the soft masses of the first two generations are disunihed from the 
respective ones, and , of the third generation. 
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experimental sources, these parameters need to be scanned numerically. We first carried out 
scans for the NNUHM, using the MultiNest-v2.18 [56] package, which is linked to the public 
code NMSSMTools-v 4.2.1 [57-60] to calculate the SUSY mass spectra and branching ratios 
(BRs) as well as the Higgs boson signal rates. The signal rate is defined, for the X decay 
channel of a given NMSSM Higgs boson. Hi, as 

nX _ ^{99 Hj) X BK{Hi —)• X) 

cr(fl'5 —^ ^Sm) X BR(/isM —> X) ’ 

where /ism is the SM Higgs boson with the same mass as Hi. 

We consider two distinct scenarios for the Higgs sector: one scenario where the Hobs 
corresponds to the lightest NMSSM Higgs boson. Hi, and one where the H^hs corresponds 
to the heavier H 2 . In the latter scenario, there is thus a Higgs boson with mass less than 
125 GeV that remains to be discovered. Both of these scenarios are easily realisable in 
significant portions of the NMSSM parameter space [61-66]. Separate scans were performed 
for the NNUHM for these two scenarios. A Hi, with i = 1, 2, is identified with the H^hs 
by requiring it to have a mass lying in the 122 — 128 GeV range^ and consistent with 
the experimentally measured = a{pp —)• Hobs X)/(7{pp —)> hsu —^ X) within ilu 
error, for each X. This condition assumes that the inclusive pp cross section for Higgs boson 
production at the LHG can be approximated by the dominant gluon fusion channel. Note, 
however, that in the rare cases when a{gg —)> Hobs —^ X)/^{99 ^SM —>■ X) for a given X is 
also provided by the experimental collaboration, we use it instead as for comparing 
with. 

The most recent publicly released values of measured by the ATLAS and GMS 
collaborations are given in table 1. Instead of considering the two sets of experimental 
measurements of the same observables separately in two scans for the given Hobs case, we 
delineate an ‘optimal allowed range’ of each H^g. This range spans min[/r^ — |1 cj|^ (ATLAS), 
p^ — |lcr|^ (CMS)] to viiax[p^ + |lcrp (ATLAS), p^ + |l(j|^ (CMS)], where p^ denotes the 
central value of the measurement and |1 ct|^ and |1 (t|^ imply, respectively, the positive and 
negative error on it, as long as p^ — \\a\^ for one experiment is not higher than p^ — \ 1 (t|^ for 
the other. This range is also given in the last column of table 1 for each X. One can notice 
in the table that while p^'^ from the two collaborations are mutually consistent, the ATLAS 
measurement of p^^ is considerably higher than the CMS one. Note also that since the WW 
and ZZ decays of a given Hi depend on the same VVHi reduced coupling, NMSSMTools 
computes a unique value, of the signal rates for both these channels. Therefore, our 

defined optimal range of covers p^^ — |lo-|^ from CMS to p^^ + llal^ from ATLAS. 

We ignore the p^^ measurement here since, aside from the fact that the error on it is very 
large, it does not take into account the gluon fusion mode for the production of Hobs- 

During our scans the Yphysics observables were subject to the following constraints. 

• 2.63 X 10“^ < BR (5 ^ Xs 7 ) < 4.23 x [71-73], 

• 0.71 X 10-^ < BR(Bu ^ Tu) < 2.57 x 10“^ [71-73], 

• 1.3 X 10-9 < BR(Bs ^ p+p-) < 4.5 x 10-9 [74-76]. 

The theoretical evaluation of these observables was carried out using SuperIso-v 3.4 [77] 
and the above ranges are the ones allowed at the 95% confidence level (CL), as suggested 

^The extended mass range is to allow up to 3 GeV uncertainty in the theoretical prediction of the mass of 
the assumed Hoha- 
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A 

fi^(CMS) [67] 

(ATLAS) 

Allowed range 

Observed 
Hobs ~ Hi 

Khs range 
Hobs ~ H2 

77 

1.13 ±0.24 

1.17 ±0.27 [68] 

0.89 - 1.37 

0.91 - 1.1 

0.89 - 1.12 

ZZ 

1.0 ±0.29 

[69] 

0.71 - 1.31 

0.95 - 1.05 

0.88 - 1.05 

ww 

0.83 ±0.21 

1.09^^;^? [70] 

TT 

0.91 ±0.28 

[69] 

0.63 - 1.9 

0.9 - 1.01 

0.63 - 1.06 


Table 1: Higgs boson signal rates measured by the CMS and ATLAS collaborations, their 
ranges enforced in the scans and those observed for the good points from the scans corre¬ 
sponding to the Hobs = Hi and Hobs = H 2 cases. 


in the manual of the package. Each SUSY point was also required to be consistent with 
the LEP and LHC exclusion limits on the other, non-SM-like, Higgs bosons of the model, 
as tested using HiggsBounds-v4.1.3 [78-81]. The satisfaction of the LEP limit on xt 
mass and the perturbativity constraints on the various Higgs boson couplings are ensured 
by NMSSMTools intrinsically. While it also tests each point against the squark/gluino 
and slepton searches from the LHC, for added surety we further required the scan to allow 
a point only if it gave a gluino mass larger than 1400 GeV. Einally, we enforced an upper 
limit on the Xi thermal relic density, to be less than 0.131, thus allowing up to -|-10% 

error in its theoretical calculation (and in this way taking into account also the experimental 
error), performed by the package MigrOMEGAs-v4.1.5 [82], given the measured value of 
0.119 from the Planck telescope [83, 84]. No lower limit is imposed on in view of 

the possibility that Xi can make up 100% of the DM in the universe via, e.g., non-thermal 
production [85], despite yielding an under-abundant thermal relic density. There is also the 
alternative possibility of the Xi constituting only a fraction of a multi-component DM (see, 
e.g., [86]). 

Our scans showed that the NNUHM serves as a good approximation for the general 
NMSSM for the Hobs = H 2 scenario, since we obtained a sufficiently large number of points 
with the desired neutralino composition. We point out here that in this scenario, the H 2 
can obtain a large tree-level mass in a more ‘natural’ way, i.e., without requiring very large 
corrections from the SUSY sectors and hence forcing sparticle masses to be large, for its 
physical mass to lie near 125GeV. This happens for A > 0.1 and tan/3 not too large. The 
scanned ranges of the input parameters were chosen in light of these observations, and are 
listed in table 2(a). Also, in this scenario the Xi^ typically a higgsino-singlino mixture, is 
generally light but can still give the correct The singlino and the higgsino masses 

are pushed towards smaller values for two main reasons: a) A is generally large and b) the 
masses of singlet-like Hi and Ai, which scale with ks, are smaller than the Hobs rnass, thus 
forbidding k and s from taking larger values. 

The Hi can play the role of Hobs over the complementary, and hence comparatively 
much larger, NMSSM parameter space, where tan/3 is not small simultaneously with large A. 
However, in the NNUHM case, a nested sampling of the NNUHM parameter space (as was 
done in our scans) tends to direct the scans towards regions yielding Hi solutions which are 
very (MS)SM-like, owing to the collider constraints imposed, given that sufficiently large niQ 
and/or oq are allowed. In such regions, it is much harder to obtain a Xi with > 0.05, since 
A is typically small and also since the relic density constraint is more readily satisfied by a 
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NNUHM parameter 

Scanned range 

mo (GeV) 

500 - 2000 

mi /2 (GeV) 

300 - 1000 

ao (GeV) 

-5000 - 0 

tan j3 

1 - 25 

A 

0.15 - 0.7 

K 

0.001 - 0.4 

/^eff (GeV) 

100 - 200 

AIa (GeV) 

-1000 - 10000 

A. (GeV) 

-500 - 2500 


(a) 


NMSSM-14 parameter 

Scanned range 

Vg, (GeV) 

1000 - 10000 

Mfj, (GeV) 

1000 - 10000 

(GeV) 

1000 - 10000 

(GeV) 

2000 - 10000 

Mj(GeV) 

500 - 10000 

Ml (GeV) 

500 - 10000 

M 2 (GeV) 

500 - 10000 

Ao (GeV) 

-25000 - 0 

(GeV) 

100 - 2000 

tan/3 

1 - 70 

A 

0.001 - 0.7 

K 

0.001 - 0.7 

Aa (GeV) 

0 - 25000 

A. (GeV) 

-25000 - 0 


(b) 


Table 2: Scanned ranges of the input parameters for (a) the NNUHM with Z^obs = H 2 and 
(b) the NMSSM-14 with i^obs = Hi. 


gaugino or a higgsino or their admixture. In our scans of the NNUHM for the H^hs = Hi 
scenario, we found only a handful of points with a non-vanishing singlino component. We 
therefore carried out numerical scans of the more general NMSSM-14 for this scenario only, 
in order to obtain a larger number of points for which bad a varying but sizeable singlino 
component. The scanned ranges of the NMSSM-14 parameters are given in table 2(b). 

In figure 1(a) we show as a function of the Xi mass, m^o, for the good points 

obtained from our scan for the NNUHM with Z^obs = H 2 ■ The colour code in the figure shows 
the distribution of Zg. The maximum Vt-oh? seen for the points is 0.131, which allows for, as 
noted earlier, up to -|-10% theoretical error given the Planck measurement of 0.119. A line 
corresponding to the —10% error, at U~o/i^ = 0.107, is also shown. The points above this 

A1 

line are considered to be in agreement with the Planck measurement and are henceforth we 
refer to as the ‘Planck-consistent’ ones. We note an appreciable density of such points with 
60GeV < m^o < 100 GeV, with Zg for them increasing with decreasing rn^o. These points 
are particularly interesting since, as noted from figure 1 in [87], Xi solutions lying within this 
mass range and also giving the correct are almost unavailable in the MSSM. There 

Al 

are another two distinct sets of Planck-consistent points, one with m~o — 32 GeV and the 

Al 

other, also an NMSSM-specific one, with m^o < 10 GeV. Zg for these points is always very 
large, m^o in this scenario does not exceed ~ 120 GeV, for reasons discussed above. 

A 4 

Figure 1(b) similarly shows the distribution of m^o for the good points from our 
NMSSM-14 scan for the Ffobs = Hi scenario. Here also a very interesting m^-o range can be 

Al 

seen for the Planck-consistent points. Neutralinos giving a good and having a mass in 

Al 

the ~ 500 — 1000 GeV range are also nearly impossible to realise in the MSSM, again accord¬ 
ing to figure 1 in [87]. Similarly to the Hohs = H 2 scenario above, for the Planck-consistent 
points the singlino fraction increases as mr-o drops, even though the typical value of the latter 

Al 

is much larger. The large density of points with m^o > 1 TeV and a fairly small Zg generally 
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Figure 1: Relic density of the x? as a function of its mass for (a) the R^obs = H 2 case and 
(b) the ifobs = Hi case. The colour code corresponds to the singlino fraction. 


is owing to the fact that a highly higgsino-dominated neutralino can easily satisfy the 
constraint for such masses. In fact none of the points with m^o > 1.4 TeV had Zg larger than 
0.05, which is also evident from the figure. 

Overall, one notices in the two figures the trend that, for a given a larger singlino 
fraction results in larger VL-oh?. This is expected, since a Xi with a small Zg (implying 
a large higgsino fraction) has a higher interaction strength and hence a higher annihilation 
rate compared to one with the same mass but a comparatively larger singlino fraction. There 
are thus some discrete Xi mass ranges for which the obtained is consistent with the 

Planck observation. In between these regions no Planck-consistent neutralino solutions were 
obtained due to the fact that for the corresponding masses the non-singlino fraction becomes 
too small to obtain sufficient DM (co-)annihilation. 

In the following, we shall not distinguish between points belonging to the NNUHM or 
to the NMSSM-14, since we are only concerned with the composition and the mass of a given 
x\- Whether a y? corresponds to the Hohs = Hi scenario or to the Hobs = H 2 scenario can 
be easily inferred from its mass. 

4 Dark matter in the Sun 

The description in this section follows the notation of [82] . The number of neutralinos in the 
Sun as a function of time, N^{t), is given by an evolution equation [88] 

^^=Cx-4«JVx(P-E,iV,(t), (4.1) 

where is the neutralino capture rate (assumed to be constant), E^N^{t) is the evaporation 
rate, and is the annihilation rate of a pair of neutralinos. Here A-^y- is defined as 




















where {av) is the averaged annihilation cross section times velocity, and is the effective 

volume of DM in the Sun [82] . Evaporation can be considered negligible for DM candidates 
above a few GeV [89], in which case oc Thus, provided that the 

capture and annihilation rates are large enough, an equilibrium will be reached after a typical 
time scale of a few times ^ after which capture and annihilation balance each 

other out and the number of the DM particles is roughly constant, iV^ ~ sfC^^JA^. The 
annihilation rate is then given by the capture rate alone, T^^ = Ay.y.N‘^12 = C^f2. For the 
Sun, most models lead to equilibrium, but MicrOMEGAs does not assume this is the case 
and solves the evolution eq. (4.1) numerically. The capture rate depends on the cross section 
for WIMP-nucleus elastic scattering and on the local density and velocity distribution in the 
Galaxy [88]. The default value of the local DM density, po, in MicrOMEGAs is 0.3 GeV/cm^ 
and a Maxwellian velocity distribution is assumed. 

The annihilation of a pair of WIMPs can lead to many different final states which decay 
further and eventually yield neutrinos. The energy spectrum of the neutrinos produced in 
the centre of the Sun is computed by MicrOMEGAs using the results of [90], for each point 
in the NMSSM parameter space considered. The resulting neutrino flux at the Earth is then 
obtained by MicrOMEGAs by computing the appropriate BRs using the NMSSM model 
(which in turn calls NMSSMTOOLS). The number of events expected in IceGube is then 
calculated by convoluting this flux with the detector effective area and angular resolution, as 
explained in the following sections. 

5 The IceCube and PINGU detectors 

The IceGube neutrino detector [40] at the South Pole consists of one km^ of very clear ice 
instrumented with 5600 optical modules arranged in 86 vertical strings. The optical modules 
register the Cherenkov radiation emitted by charged leptons produced in high-energy neutrino 
interactions in the ice. For the purpose of this paper, we are particularly interested in muons 
from the interactions of muon neutrinos, because such muons are observed as tracks giving 
a good angular resolution. Although the IceCube geometry has been optimised to detect 
ultra-high energy (> TeV) neutrinos from potential cosmic sources with sub-degree angular 
resolution, its angular response of a few degrees at 0(100 GeV) is still adequate to perform 
directional searches for DM. Unfortunately the energy of muon-neutrino induced events can 
not be precisely measured since the long muon track might not be contained within the 
detector volume. Electron neutrinos and tau neutrinos, on the other hand, leave localised 
energy depositions in the detector that, while offering much worse angular resolution, allow 
for a better energy measurement since they are contained in the detector. However, since 
we are interested in neutrinos originating in the centre of the Sun, it is important to be able 
to identify the direction, and we are therefore going to consider only muon neutrinos in the 
following. In the centre of IceCube there is a more densely instrumented sub-array called 
DeepCore, which is more sensitive to lower energy neutrinos, having an energy threshold 
around 10 GeV [91]. As we will see below, it is very important to have sensitivity to low-mass 
DM. In addition to the IceCube and DeepCore detectors, we consider a proposed extension 
called PINGU, a low-energy in-fill of about 40 additional strings with closer optical module 
spacing, which will have an energy threshold for neutrinos of a few GeV [41]. 

The efficiency of a neutrino telescope to a given signal is parametrised by the effective 
area for muon neutrinos. A®®, which is a function of the neutrino energy. The effective area 
is the equivalent area where the detector can detect a neutrino with 100% efficiency, and 
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includes the neutrino-nucleon interaction probability, the energy loss of the produced muon 
from the interaction point to the detector and the detector trigger and analysis efficiency. For 
the IceCube/DeepCore analysis we use the detector effective areas presented in [92], which 
were obtained from detailed simulations of the detector response within the framework of a 
search of IceCube data for DM in the Sun [93] . Below 40 GeV we use the PINGU effective area 
presented in [94]. We then have three practically independent detectors, together covering a 
large neutrino energy range, with little overlap between their individual coverages, as shown 
in figure 2 below. 


6 Signal and background estimations 


6.1 Expected signal from the Sun 

We want to calculate the number of events that IceCube, DeepCore and PINGU can expect 
for the different NMSSM model points we analyse, and compare it with the number of 
expected background events from the direction of the Sun. The neutrino flux <I>, obtained 
using MicrOMEGAs from the NMSSM model, is parametrised as 


4 > = 


dA dt did ’ 


( 6 . 1 ) 


with neutrino energy Ei,^, area on the Earth A, time t and solid angle il of the sky. D 
is defined in terms of the polar angle 6 and the azimuthal angle ip as dil = dipdO sin0 . 
MicrOMEGAs returns the neutrino and anti-neutrino fluxes separately. Even if we take 
“neutrino flux” to denote both neutrino and anti-neutrino flux in what follows, we treat 
neutrinos and anti-neutrinos separately in our calculations and use the corresponding cross 
sections when needed. Even though the signal from DM annihilation is independent of t, 
the background depends on the position of the Sun and, therefore, on t. For simplicity, we 
choose the coordinate system with 0 = 0° in the direction of the Sun. For a given detector, 
the total number of events observed from an incident neutrino flux is given by 


/ roc |^2^T rOcut 

dtj^ dE^^J^ dp>j^ deAueAf^{E,^)^{E,^,e,p,,t), ( 6 . 2 ) 

where A®® is the effective area. The integral over t is done over the live time of the detector, 
as will be specified later. In our case of a search for DM annihilation in the Sun, the angular 
integral is performed over a cone that points to the Sun’s position with a cut on the maximum 
polar angle Ocut to take into account the angular resolution of the detector. at the detector 
location (Earth) will be modified by nentrino energy losses as well as oscillations on their way 
from the Sun. MicrOMEGAs takes care of these effects, providing the expected neutrino 
flux at the Earth for each signal model point considered. In a similar way, events from 
anti-neutrinos are calculated. 


6.2 Atmospheric background 

Neutrinos created in the atmosphere by cosmic ray interactions (see e.g. [95-98]) constitute an 
irreducible background to the signal, all year round. On top of this background, atmospheric 
muons produced in the same cosmic ray interactions provide an additional background that 
can be dealt with differently depending on the time of the year. During the antarctic summer 
season, when the Sun is above the horizon, neutrinos from the Sun are downgoing. Any search 
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for a solar signal is then plagued by the overwhelming atmospheric muon background, which is 
about seven orders of magnitude larger than the atmospheric neutrino background, depending 
on the declination [96, 98]. It is only through vetoing techniques (considering events starting 
inside the detector volume, which must be caused by a neutrino) that IceCube has managed 
to perform searches for DM in the Sun in the austral summer also. The price to pay is 
a significantly reduced detector volume, since part of the outer strings of the detector are 
used for tagging incoming atmospheric muons. Such an analysis requires access to detailed 
detector event information, which is beyond the scope of this work. We will thus concentrate 
on a “winter”-type analysis in the rest of the paper by using the corresponding effective areas 
from [92]. 

During the antarctic winter season, the Sun is below the horizon at the South Pole and 
neutrinos coming from it are upgoing, meaning that they pass through the Earth before arriv¬ 
ing at IceCube. This provides a hlter against the atmospheric muon background. The only 
physical background then is the atmospheric neutrino flux. There is, however, an additional 
background component consistent of downgoing atmospheric muons misreconstructed as up- 
going. This is a background that can be minimised by clever data analyses, but a residual 
contamination of a few percent can always be present in the final analysis. We will assume 
in the rest of the paper that such misreconstructed atmospheric muon contamination can be 
reduced to an insignificant level in a real event-by-event analysis, and we will not consider 
such a background further in our calculations. Our results are then slightly optimistic, but 
we note that a few percent increase, e, in the background would worsen the significances 
presented in section 7 just at the level of 1/ y/l + e. 

For our estimations of the PINGU sensitivity, on the other hand, we will not make any 
distinction between winter and summer periods. This is because we assume that using the 
surrounding IceCube detector as a veto, downgoing atmospheric muons can be effectively 
rejected in the summer period and the detector can be used all year round in the same way. 

We calculate the expected background from state-of-the-art atmospheric neutrino flux 
calculations. At the energies relevant for our analysis, this background is dominated by the 
conventional neutrino flux from decays of pions and kaons. The prompt flux from charmed 
hadrons is only relevant for much higher energies. The conventional flux depends on the 
neutrino energy and the zenith angle. The flux prediction by Honda et al. [98] has been 
modified by Gaisser [99] to account for the cosmic ray knee and agrees well with the measured 
flux. It is used by the IceCube collaboration as a theoretical standard prediction for the 
conventional flux. We use this flux for our background estimates through the implementation 
in the code NeutrinoFlux [100]. The calculation proceeds along the same lines as for the 
signal estimation, that is, through the convolution of the atmospheric neutrino flux with the 
effective area/volume of the detector as in eq. (6.2). 

7 Results and discussion 

The predicted spectra from DM annihilations in the Sun are obtained from the NMSSM 
model points via MicrOMEGAs, and the expected signal is calculated using the effective 
areas of the detectors as indicated in eq. (6.2). We consider the three independent detectors 
IceCube, DeepCore and PINGU, as noted above and illustrated in figure 2, each covering a 
different range in neutrino energy. Figure 2 also contains examples of the spectra resulting 
from DM annihilation in the Sun. A few spectra with the largest signal were selected to 
illustrate the general features of the distributions. The irregular structure of some of the 
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Figure 2: Examples of the spectra (left axis) from some selected model points which 
result in a large statistical significance. The non-smooth features in some of the spectra 
originate from DM annihilation into WW or ZZ. Also shown is the energy dependence of 
the PINGU (orange), IceCube (blue) and DeepCore (violet) effective areas (right axis). 


spectra around 30 GeV is due to the annihilation channel x° x° —)• W~^W~, which can result 
in such characteristic features if the DM is heavy enough so that this channel is possible, but 
light enough for the spectrum to not get smoothed out [90]. 

For a predicted number of signal events Ng from a model and a given number of back¬ 
ground events we estimate the significance S by the formula [101] 

where the last approximation is valid for Ng <C N^,. The complete expression above yields a 
lower S compared to the approximate one in the case of large Ng and the difference between 
the two amounts in our case to a maximum correction of 20%. 

We integrate the signal and background over the actual live times of the detectors during 
one calendar year. As pointed out earlier, PINGU is assumed to take data during the whole 
year, since we assume that IceGube can be used as an efficient veto. IceCube and DeepCore 
are assumed to take data only during 152 days of the winter time as specified in table 7.1 
of [92]. This lets us compare the performance of the different detectors in a more fair way. 

We cut on the angle 6 between the reconstructed track and the direction of the Sun as 
outlined in eq. (6.2) to take into account the reconstruction error for the neutrinos from the 
point-like source. For the signal, we assume that 9 is normal-distributed around 0° with a 
median given by figure 7.17 in [92]. S is then maximised over ^cutj with the typical optimum 
values being 3° for IceCube and 8° for DeepCore, while we use 20° for PINGU. The exact 
optimum value depends also on the specific NMSSM model point, since the angular resolution 
depends on the neutrino energy. 

In the context of the current direct and indirect DM searches, figure 3 shows the spin- 
independent and spin-dependent DM-proton elastic scattering cross sections, denoted by 
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Figure 3: Scatter plots of the model points showing CTp^si (left) and cTp^sD (right) in terms 
of the Xi mass, with the colour code corresponding to The top panels contain all 

the good points from the scans, while the lower panels show only the model points with the 
predicted lying within ±10% of the Planck measurement. 
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(Tp^si and (Tp^sD, respectively, as functions of rn-o, with the colour-code indicating the value 
of Q.^oh? predicted for each model point. The top panels show the complete set of good 
points obtained from all our scans, while the lower panels contain only the Planck-consistent 
(i.e., with 0.107 < < 0.131) subset of these points. The experimental limits shown in 

the left panels correspond to the current 90% CL exclusion limits from XENONIOO [30] as 
well as the projected 90% CL limits for XENONIT [102]. In the right panels we similarly 
show the 90% CL exclusion limits from SIMPLE [103] and IceCube [93], as well as the 
projected 90% CL exclusion limits for PINGU [41]. The current limits from IceCube are 
shown as a guideline, but they should be interpreted with care. They were obtained under 
some simplifying assumptions - 100% annihilation rate into either bb (soft) or W^W~ (hard) 
- and therefore can not be taken to exclude all model points lying above the corresponding 
lines. A real NMSSM model point has a mixture of annihilation channels and these limits 
would be relevant for a parameter point only if it predicts a neutralino annihilation rate 
almost exclusively into either bh or W~^W~. 

As can be seen from this figure, a fraction of the considered parameter points with 
m^o ~ 100 GeV and m^o ~ ITeV is already excluded by XENONIOO, and a large fraction 
of the remaining such points can be tested with the upcoming XENONIT experiment. The 
XENONIT and IceCube exclusions are, however, complementary, which we will come back 
to below. A note of caution here is the well-known fact that the direct and indirect detection 
limits on SUSY models are subject to different uncertainties from particle and nuclear physics 
and astrophysics. In particular, the direct detection limits are affected by large uncertainties 
from nuclear physics [104-106] (for recent analyses of the impact of such uncertainties on the 
cMSSM and MSSM parameter spaces see, e.g., [107-109]). Indirect searches for DM in the 
Sun are affected by several astrophysical uncertainties, mainly on the value of the local DM 
density, the assumed velocity distribution and velocity dispersion, although they have a more 
moderate effect on the final limits [110, 111]. Therefore, while we consider the points that 
lie above the line corresponding to XENONIOO in figure 3 to be excluded, we shall anyhow 
include them in our discussion below, but separately from the allowed points, in order to 
illustrate the full complementarity of the direct and indirect searches. 

The impact of IceCube, DeepCore and PINGU on the higgsino-singlino DM predictions 
in the NMSSM is illustrated in figures 4 and 5. The figures show the significance S defined 
in eq. (7.1) as a function of the Xi mass for the NMSSM model points after 1 calendar year. 
Also drawn in the figures are the lines corresponding to discovery (5a} as well as to exclusion 
(2a} after 1 or 10 calendar years, the latter aiming to highlight the model points that have 
not been tested by any of the detectors yet but will be after 10 calendar years of data-taking. 
Note that the models above the 5a line in the plots would have probably been seen in the 
IceCube/DeepCore analysis presented in [93] which, even when optimised for the MSSM Xij 
does not show any other specific model dependencies, and is general enough to be sensitive 
to the NMSSM x)* as well. Such points are thus in principle already ruled out. The colour 
code in both the figures shows Z^, and the model points have been sorted in such a way that 

a point with larger Zg lies on top of another one with the same m^o and S. 

•^1 

Figure 4 shows the best significance, Ssest, obtained for a given model point by either 
IceCube or DeepCore. Results for PINGU will be shown later in a separate figure because 
PINGU has not yet started to collect data. The left panels contain only the Planck-consistent 
points while the right panels contain the Planck-inconsistent ones, i.e, for which VL~oh? < 
0.107. The top panels show points that are not excluded by XENONIOO, while the lower 
panels contain the complementary set, i.e., models disfavoured by XENONIOO. One sees in 
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Figure 4: Overview of the maximum significance for each model point obtainable from 
the IceCube or DeepCore detector. The cut 9cut on the angle between the tracks and the 
direction of the Sun is optimised for each model point. The values for the significance 
represent 1 calendar year of data-taking. IceCube and DeepCore are assumed to operate 
only during 152 days of winter, as described in the text, but PINGU, shown in later figures, 
operates independent of season. The colour code indicates the singlino fraction. 


the figures a high sensitivity to a large number of points with m^o ~ 100 GeV, which can 
partly be traced back to the DM annihilation channel XiXi W~^W~ yielding a harder 
spectrum for such a Xi mass [90]. In fact, a large fraction of the m^o ~ 100 GeV points 
in each panel lie above the 2a line shown. The top panels are of particular interest in this 
regard, since they indicate that many of these points, potentially ruled out by IceCube after 
just one calendar year of data-taking, were not probed by the XENONIOO experiment. 

Figure 5, on the other hand, shows the performance of each detector individually, but 
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Figure 5: Expected significances at the IceCube (left), DeepCore (middle) and PINGU 
(right) detectors. We note that DeepCore already rules out a large fraction of the parameter 
space points with m^o ~ 100 GeV with one calendar year of data-taking (cf. hgure 4), and 
PINGU has the potential to reach lower-mass points than those accessible with DeepCore 
within a few years of runtime. All points shown pass the XENONIOO limits. 


for model points that satisfy only the upper limit on and are allowed by XENONIOO. 

We see that it is DeepCore (middle) and PINGU (right) which contribute more strongly to 
rejecting a large number of the m^o ~ 100 GeV points. As can be seen in figures 2 and 3, 
the XiXi W^W~ annihilation channel is well covered by the DeepCore detector for this 
m^o range. Almost all the model points belonging in the 200 GeV < m^o < 1400 GeV strip 
noted in the left panel (for IceCube) of figure 1(b) he just below the line corresponding to 
2a significance after 10 years and will thus not be accessible within this timescale. Note 
that the inclusion of PINGU in our analysis is mainly to estimate its sensitivity to the 
model points with m~o < 20 GeV. The results show that PINGU is indeed sensitive to 

A1 

Xi in the 0(10 [GeU]) mass range, although it needs a few years of accumulated livetime. 
Dedicated event-based likelihood analyses, on the lines of the one presented in [112], are 
known to improve the sensitivity of an experiment since they use more information in a more 
comprehensive way, and can enhance the prospects for reaching such points with higher 
significance in a shorter timescale. To evaluate the effect of such an approach is, however, 
beyond the scope of this paper. 

In figure 6 we show again the S'sestj but with the singlino fraction now on the x-axis 
and the DM mass given by the colour code. The figure further illustrates that the IceGube 
experiment can have a reasonable sensitivity to an abundance of points with a wide range of 
singlino fractions. 

Finally, in figure 7 we highlight the complementarity of DeepCore/PINGU and the 
XENON direct detection experiment. The figure shows the m~o ~ 100 GeV model points 
which can not be excluded by the projected XENONIT 90% GL limits in a light shade on 
top of the ones that can be excluded, shown in darker points. As seen earlier, a subset of 
the points is already excluded by the available DeepCore one-year data. We note that more 
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Figure 6: Best significance (from either of IceCube, DeepCore or PINGU) obtained for 
model points allowed (left) or excluded (right) by XENONIOO, against the corresponding Zg 
and m^o. 


parameter space points, some of which will not be probed by XENONIT, will be testable with 
each subsequent year of data-taking by the DeepCore detector. The PINGU detector, when 
it enters operation, will also be able to test model points reachable neither by XENONIT 
nor by DeepCore due to the low m~o, thus illustrating the complementarity of PINGU and 

A1 

DeepCore. 

8 Summary and conclusions 

In this article we have analysed in depth the potential of the IceCube/DeepCore detectors 
and the proposed low energy extension, PINGU, to probe the higgsino-singlino sector of the 
NMSSM. We have presented our findings from the NMSSM parameter space scans, which 
concentrated on model points predicting a lightest neutralino with a non-vanishing singlino 
fraction, and hence a non-MSSM-like dark matter candidate. In these scans we required such 
points to survive the most important experimental constraints from Higgs boson searches as 
well as from 5-physics and relic density measurements. 

We have used the spectra from DM annihilation in the Sun to test the model points 
of our interest. We have evaluated the signal on a point-by-point basis in the NMSSM and 
taken into account the expected atmospheric background in the detectors considered. We 
have then estimated the statistical significance that can be obtained for each model point by 
each detector and found that a subset of all the points obtained from our scans have already 
been ruled out by the one-year data from the IceGube experiment, while more points will be 
accessible with each subsequent year of data-taking. We have found that PINGU will indeed 
have sensitivity to neutralino masses of the order of 10 GeV, although the detector should 
accumulate several years of livetime to reach a significance of 2ct for some of the points with 
a low mass. 
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Figure 7: Number of calendar years required to reach S' = 5cr as a function of the Xi 
mass. The left plot shows only points where PINGU gives the best significance per calendar 
year, whereas the right plot illustrates in a similar way the case of DeepCore. Points not 
reachable by XENONIT with the projected sensitivity are shown in a lighter shade on top of 
the remaining points. We conclude that DeepCore and PINCU complement each other, and 
can access model points that are out of reach of the projected XENONIT sensitivity. The 
number of years has been extrapolated from the averaged 1-year significance. For PINCU 
we assume continuous data-taking over the calendar year. 


We have also emphasised the complementarity of IceCube/DeepCore/PINCU searches 
to the ton-scale direct detection experiments, by highlighting parameter space regions that 
may not be probed by the latter but will be accessible at the former. 
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